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Abstract 



A practical criterion for the irreducibility (with respect to integration by part iden- 
tities) of a particular Feynman integral to a given set of integrals is presented. The 
irreducibility is shown to be related to the existence of stable (with zero gradient) 
£T) | points of a specially constructed polynomial. 
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! The growing accuracy of high energy physics experiments demands the calcu- 

lation of higher order quantum corrections. The latter, in turn, are expressed 
through multi-loop Feynman integrals. Most of the methods for their calcu- 
lations are designed for the integrals with a monomial of momenta in the 



Ph. numerator. On the other hand, in practice one can face thousands and even 

millions of such integrals in one particular physical problem. Fortunately, it is 
possible to reduce this number to around a few dozens by using integration by 
part identities [1]. These identities can be solved in some cases by construct- 
ing recursive algorithms [1,2] or by "brute computer force" with the Laporta 
algorithm [3]. In any case, it is useful to know if the given integral can be re- 
lated to simpler integrals (with some denominators missing). A corresponding 
criterion of irreducibility was proposed in [4]: the irreducibility was related 
to the existence of a special solution of the recurrence relations; a recipe to 
construct such solutions was suggested. Unfortunately, the recipe in the form 
of [4] required some "hand work" which becomes very undesirable for the real- 
life problems (too many cases to be considered). So in this paper we propose 
a more practical version of this criterion which amounts to the study of zero- 
gradient points of some polynomial and hence can be easily implemented with 
a computer. We start with examples, then formulate the criterion and finally 
present the list of the irreducible four-loop massless propagator integrals found 
with the help of our criterion. 
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1 Examples. 



Consider two massless 4-loop propagator integrals: 




"cube" "ladder" 

The "ladder" integral consist of triangles and thus, according to "triangle rule" 
[1], can be related to the combination of simpler integrals; the reducibility of 
the "cube" diagram is unclear. To study this question we are going to associate 
with every multi-loop integral some polynomial, whose properties reflect its 
irreducibility. 

Let us calculate the determinant P = det(piPk) {pi = q,m,p,k,l) under the 
condition that all lines are "on mass-shell", that is corresponding squared mo- 
mentum are equal to zero. In total there are 15 scalar products ("external" q 2 
and 14 "internal" involving loop momenta) with 11 constraints, so we can ex- 
press 11 "internal" scalar products through others. For "cube" it is convenient 
to choose pq, Ip, Iq, q 2 as independent: 

Pcub = -(2 Ip ■ lq + 2 Ip ■ pq + Ip ■ q 2 - 2 Iq ■ pq ) ■ Ip ■ Iq ■ pq. 
For "ladder" let us choose kl,kp,lp,q 2 as independent: 
Piad = (kl -kp-lp- q 2 /2) ■ kl ■ kp. 

For each polynomial let us check if there are points where all derivatives with 
respect to the "internal" scalar products are zero, but the polynomial itself 
and the determinant of second derivatives are not zero (in short "non-zero 
stable points"). The "ladder" polynomial has no such points, and the "cube" 
has one point: Ip = —9q 2 /20,pq = —3q 2 /10,lq = — 3g 2 /10. As we will see 
later, this means irreducibility of the "cube" integral. 



2 Criterion in general form. 

Now let us explain why non-zero stable points of some polynomial are related 
to the irreducibility. Assume that there is a set of equations which relate the 
function B(n) (keeping in mind Feynman integrals as functions of degrees of 
propagators and irreducible numerators) with arguments shifted by ±1: 

R(r,I + )B(n) = 0, (1) 
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where (n) (and others underlined arguments) is (ni, .., njv), I~B(..,n a , ..) = 
B(..,n a — 1, ..) and 1+ £?(.., n a , ..) = n a B(.., n a + 1, ..). Since I~ and I + do not 
commute ([I~, = <5 a b), for defmiteness we place all I + to the right of I . 

We want to check the possibility to represent, using the relations (1), a given 
B(n^) as linear combination of {B(n^), i e (1, k)}: 

k 

0lB(nW)-ZfiB(n®). (2) 
i=i 



The relation (2), if it exists, should be a combination of the relations (1) for 
various values of n, so any solution of (1) should also obey (2). Suppose we 
are able to construct the special solution of (1) s(n) with the properties 

s(n (0) )^0, s(n«) = for i e (1, k). (3) 

This s(n) evidently cannot obey (2) and hence (2) cannot be the combination 
of (1). So we got a sufficient criterion of irreducibility in general form [4]: the 
recurrence relations (1) cannot relate the B(n^) to linear combination of the 
set {-B(n^), % G (1, k)} if a partial solution of (1) with properties (3) exits. 



3 Formal solutions of the recurrence relations. 

Let us construct the solutions of the relations (1) as formal series in 1/D: 

k=0 

To be precise, we will present an algorithm for calculation of the B^ such 
that for any n and arbitrary large N, sufficiently long partial sum of (4) will 
obey (1) up to terms 1/D N . 

As an example, let us consider the recurrence relations 

= (k - 2D + 2)B k+l + 2a(k -D + l)B k + kbB k _ t (5) 

generated by the integral 

oo 

B k = [ x k (x 2 + 2ax + b)~ D dx, k > 0, k + 1 - 2D < 0, b > a 2 . (6) 



3 



After the substitution x — xj y/D — a integral (6) can be represented as 



oo _2 

B k = f x k (%^ + b-a 2 r D dx. 



D 



-oo 



The B k obeys the recurrence relations 

= (k(b - a 2 )B k ^ - 2B k+1 ) + + 2)B k+1 = R° k (B) + ^R\{B). (8) 



Inserting B k = Y4I0 /D l we obtain the following equations for B®: 



= R° k (B (0) ), = R° k (BV +1) ) + Rl(B {l) ) I > 0. 

If b — a 2 7^ 0, these equations have nontrivial solution and we obtain the 
desired solution of the (5) as formal series in 1/D. Note that we can forget 
about integral (6) and study the relation (5) with the help of substitution 
(7). In particular if b < a 2 the integrand of (6) has poles on the real axes, 
but the substitution (7) does not depend on the definition (and even on the 
convergence) of the integral. From other side, the "integral representation" 
language is more convenient, so in general case we will construct solutions 
(3) presenting necessary manipulations in this way, keeping in mind that for 
formal proofs we can translate all steps to "recurrence relations" language. 



4 Formal solutions in general case. 

To construct the integral representation for s(n) (3) let us try [5]: 

4«) = / (9) 

The action of R(I~,I + ) on (9) leads to 



R(I ,I + )s(n)= f dxi...dx N g(x)R(x,-dx) ni 

J X^ — X 



n N 
N 



I ^nl'"^ R'U, dx)g(x) + (surface terms), (10) 

J Xi ---Xj^ 
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where we use integration by part in second equation and denote by R' the 
combination of x and dx acting on g(x). Now, if one chooses g(x) so that 

&(x,d s )g(x) = (11) 



and chooses the integration contours such that the surface terms vanish, one 
arrives to the function s(n) such that R(I + , I~)s(n) = 0. 

As shown in [5], the (11) can be solved for the general case of a multi-loop 
Feynman integral with arbitrary number of legs and with arbitrary masses, and 
the corresponding g(x) can be represented as a product of two polynomials in 
x, each polynomial raised to non-integer degree. Then, the recurrence relations 
for L-loop i?-leg integrals is a partial case of the recurrence relations for 
R + L — 1 tadpole integrals [6], thus we limit ourselves by the tadpole case. 

Consider an L-loop tadpole Feynman integrals in dimensional regularization: 

= f d D Pl . . . d D p L 

where pi (i — 1, . . . , L) are loop momenta and 

D a = Y.Y.A^p l .p J -ml (a = l,...,N = L{L + l)/2) (12) 
i=i j=i 



are propagators (n a > 0) or irreducible numerators (n a < 0). Consider some 
integral B{n^) with the set of positive indexes ("lines") S: a e S if > 0. 
Normally, the integrals with less number of lines are easier to calculate, so 
let us check if it is possible to relate this integral to integrals with some lines 
missing. For that we should try to construct s(n) such that s(n) = if at least 
one of the n a( z S < 0. Eq. (9) in this case will read [5] 

= / $^n^ L - 1)/2 , (13) 



where P(x) = detrki)(^2a=i^a l ( x a + m a)) is polynomial in x and A^ 1 is the 
matrix inverse of A* 1 , i.e. J2a=iA^Af = (S ik Sji + 5u5j k )/2. In other words, 
to construct P one should calculate the determinant of scalar products of the 
loop momenta and then substitute the scalar products by the combinations of 
propagators and numerators obtained by inverting (12). 

Let us choose the integration contours for x a€ s as small circles around zero. 
In this case, according to Cauchy's theorem, the integrations will lead to the 
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(n a — l) th coefficient in the Taylor expansion of the integrand: 



s{n) oc / [] 



x'; 



n 



0: 



ma-l 



1)! 



P(x) 



(D-L-l)/2 



(14) 



In the case when the number of lines is equal to the number of integrations 
in (13) the (14) has no integrations and (if P(2L)\ X s=0 ^ 0) the construction 
of the s(n) is completed; in the general case, we need a way to define the 
integration over the remaining x a £s- As the result of the Taylor expansion, 
the original s(n) will be linear combination of the terms of the following type: 



k a ^ s a£S 

P(x) = P(x)\ Xa£S=0 , D' = D - L - 1 - £(n a - 1). 



(15) 



Let us expand (15) in the formal series in the limit D' — > — oo. Assume there 
exists a point x° with: 

P(x°) ± 0, d ats P{£) = 0, &et(d ats d HS P(£)) ^ 0. (16) 

ab 



After the substitution x = x° + x[j \J —D' 



P{x)^ « (p(*°) - £ ^d a d b P(x°) + 0((-D>)- 3/2 )) 



a,b<£S 



2D 



P(x )^(e"^ E ^ s ^ aa ^ ( -° ) + O((- J D')" 1/2 ))- (17) 



The Gaussian integrations can be performed by the formula 

r M ir M / 2 M rl 

J ll{dx a x a )e - det (C)V2 l^dtj 6 t a =o- ^ 



Finally, expansion in D' — > — oo can be easily recalculated to D — > — oo. 

Eq. (16) defines the practical version of the criterion: if (16) is valid, then the 
eqs. (14,15,17,18) define the algorithm for calculating s(n) with properties (3) 
thus proving the irreducibility of the integral B(n^). Note that P in (13) is 
det(piPk) where scalar products are expressed through linear combinations of 
x, and hence P in (15) is equivalent to det(piPk) with lines "on-shell", as it 
was considered in the examples. 
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The equation (16) and hence (15) can have several solutions s^ p \n). This 
corresponds to possibility to find several integrals with the same set of lines 
(differ by number of dots or by numerators) which are irreducible as to each 
other, as to simpler integrals. To fix these irreducible integrals, let us choose 
some set of points in "n" -space such that the matrix = s^(n^) has 
inverse matrix M q . Then solutions M q p s^ will be "diagonal" in the sense (3) 
and hence the set will define irreducible integrals. 

The application of the criterion (16) to the four-loop massless propagator 
integrals (or equivalently five-loop tadpole with one massive line) results in 
the following list of integrals which cannot be reduced to integrals with a 
smaller number of lines (for each of the topologies depicted here only a single 
irreducible integral exists): 
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5 Final remarks. 



The criterion (16) is easy to implement in the case when the polynomial P in 
(15) depends only on a few variables (large number of lines in the integral). 
In the most easy case when number of lines is equal to number of " internal" 
scalar products the irreducibility will follow from P(0) 7^ 0. In the case of a 
small number of lines (many variables in P) it may be difficult to solve (16) 
explicitly. In this case it is instructive to check (using Groebner basis) if the 
polynomial P is reducible to polynomials <9 a ^sP. If so, then d a ^sP{x°) = 
leads to P(x°) = and there no non-zero stable points. On the other hand, 
reducibility of the integrals with a small number of lines in many cases can be 
checked directly (with the triangle rule or by checking subdiagrams). 

And finally, let us address the question what happens if (16) is not satisfied. 
Formally this case is unclear, but in practice (at least for massless propagator 
integrals up to 4-loop and for zero-scale 2-loop propagator integrals [7]) it 
corresponds to reducible integrals. 

The author are grateful to K.G.Chetyrkin, R.Harlander, J.H.Kuehn, Y.Schroder 
and O.V.Tarasov for the careful reading of the manuscripts and useful discus- 
sions. 
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